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Hypothesis: Underground hydrogen (H2) storage is a potentially viable solution for large-scale cyclic Hz 
storage; however, the behavior of Hz at subsurface pressure and temperature conditions is poorly known. 
This work investigates if the pore-scale displacement processes in H2-brine systems in a porous sand- 
stone can be sufficiently well defined to enable effective and economic storage operations. In particular, 
this study investigates trapping, dissolution, and wettability of H2-brine systems at the pore-scale, at con- 
ditions that are realistic for subsurface H3 storage. 

Experiments: We have performed in situ X-ray imaging during a flow experiment to investigate pore-scale 
processes during H3 injection and displacement in a brine saturated Bentheimer sandstone sample at 
temperature and pressure conditions representative of underground reservoirs. Two injection schemes 
were followed for imbibition: displacement of Hə with H2-equilibrated brine and with non-H;- 
equilibrated brine. The results from the two cycles were compared with each other. 

Findings: The sandstone was found to be wetting to the brine and non-wetting to H3 after both displace- 
ment cycles, with average contact angles of 54° and 53° for H2-equilibrated and non-H2-equilibrated 
brine respectively. We also found a higher recovery of H (43.1%) when displaced with non-H>- 
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equilibrated brine compared to that of H2-equilibrated brine (31.6%), indicating potential dissolution of 
H2 in the unequilibrated imbibing brine at reservoir conditions. Our results suggest that underground 
Hə storage may indeed be a suitable strategy for energy storage, but considerable further research is 
needed to fully comprehend the pore-scale interactions at reservoir conditions. 

© 2022 The Author(s). Published by Elsevier Inc. This is an open access article under the CC BY license 


(http://creativecommons.org/licenses/by/4.0/). 


1. Introduction 


The Intergovernmental Panel on Climate Change (IPCC) [1] has 
warned that the world is set to reach the 1.5 °C warming level 
within the next two decades, and that only the most drastic and 
immediate cuts in carbon emissions will prevent an environmental 
disaster. Clean energy alternatives such as H2 have the potential to 
support near-zero greenhouse gas emissions [2-7]. While there 
has been an unprecedented economic and political momentum 
towards the use of green Ho, the large volumes needed for future 
demand would require new storage facilities [8-12]. 

Underground H storage, for example in porous rocks or salt 
caverns [13-16], can provide an effective solution for short- and 
long-term storage to meet the fluctuations in energy demand. Salt 
caverns have already been used for storing natural gas and Hp, for 
example, in Cheshire, Teesside, and Yorkshire [17] in the UK; how- 
ever, they have a limited storage capacity [13] (170 TWh for all of 
Europe) [16]. An underground H; storage capacity between 250 
and 1000 TWh is required for Europe by 2050, to support the fore- 
casted demand of 2500 TWh [16]. To scale up the storage capacity, 
porous geological formations located deep underground, for exam- 
ple, containing depleted gas reservoirs [16,18,19] and saline aqui- 
fers, have the potential to provide the required storage volumes for 
H2 (750 TWh capacity available in Europe) [16]. The feasibility of 
storing, and subsequently withdrawing H3 repeatedly from such 
porous underground reservoirs to meet fluctuations in energy 
demand is currently being investigated in many countries 
[10,12,16,17,20]. 

Storage in subsurface reservoirs is regarded as economically 
viable, requiring substantially less capital cost and a lower leve- 
lized cost of Hz storage (1.23-1.26 $/kg) compared to salt caverns 
(1.61 $/kg) [21]. However, there are significant uncertainties in 
storing H, in porous geological formations [22,23], which are due 
to the limited knowledge to predict and mitigate the negative 
effects of: (i) the uncontrolled migration and unwanted trapping 
of Hz in reservoir rocks [23-25], (ii) H2 reactivity with water and 
minerals [26-28], (iii) Hə diffusion through the reservoir and 
caprock [8,29], and (iv) formation of biofilms [30,31] and their 
impacts on pore blockage and storage efficiency. All these are 
major challenges associated with underground H, storage and 
need to be studied independently in detail to obtain quantitative 
information of the contribution of each process on the effective- 
ness of H; storage in the subsurface. 

Similar to other fluid-fluid displacement processes such as CO2- 
brine and CH,-brine, capillary trapping and wettability are of 
importance to H2-brine systems [32], as the amount of fluid trap- 
ping directly affects recovery efficiency [33-35]. The cyclic nature 
of an underground H; storage system renders these phenomena 
(capillary trapping and wettability) even more critical due to the 
repeated occurrence of drainage and imbibition cycles. 

Capillary trapping mechanisms have been extensively studied 
for CO», No, and oil in various porous media [36-43]. These studies 
generally focus on the pore-scale, where the capillary trapping can 
be observed and quantified directly. However, it is not yet known if 
Hə behaves similarly in subsurface rocks, in which capillary, vis- 
cous, and buoyancy forces have the potential to play a key role 
along with other geochemical and biological processes. Most of 
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the available studies on H3 storage in rocks are focused on 
reservoir-scale numerical simulations [15,44—47], which do not 
yet provide pore-scale information. 

In addition, there is still a lack of consensus on the wettability 
state for subsurface Hz storage systems. First experimental studies 
have reported a wide range of contact angle values for Hz-brine 
systems: using the tilted plate technique, contact angles of 0° to 
50° for pure quartz and 50° to 70° for a quartz sample aged with - 
stearic acid were measured for a pressure range of 1 to 250 bar and 
a temperature range of 20 to 70 °C [35]. A clear correlation 
between increasing contact angle and increasing pressure and 
temperature was observed. It was also reported [26] that aging 
rock samples in stearic acid increased the contact angles for H2- 
brine system from strongly water-wet to intermediate-wet for 
rocks containing calcite, dolomite and shale and, from strongly to 
weakly water-wet for rocks containing gypsum, quartz and basalt. 
The measurements were performed using the captive-bubble 
method for temperature ranges of 20 to 80 °C and pressure ranges 
of 10 to 100 bar. These studies conclude that under subsurface con- 
ditions, Hz wettability is prone to be intermediate-wet for storage 
in depleted oil reservoirs. 

In contrast, contact angle values between 21.1° and 43° have 
been reported for experiments performed on Berea and Ben- 
theimer sandstone samples using the captive bubble method for 
temperature ranges of 20 to 50 °C and pressure ranges of 20 to 
100 bar [24]. The contact angle was not affected by the changes 
in temperature, pressure, salinity, and rock type; however, it was 
affected by the bubble size. Similarly, advancing contact angles 
ranging between 13° to 39° and receding contact angles between 
6° to 23° have been reported for a Hz-water system using a 
microfluidic device with varying channel widths at a temperature 
of 20 °C and pressure of 10 bar [48]. Moreover, contact angles of 
21.56° and 34.9° at two different experimental conditions of 
50 bar and 20 °C, and 100 bar and 45 °C, respectively were reported 
based on core flood experiments in Voges sandstone [49]. These 
studies indicate strong water-wet conditions under all the tested 
parameters. However, to date no study has measured in situ con- 
tact angles of Hz-brine systems at subsurface pressure and temper- 
ature conditions in a real rock. 

In situ measurements have been made possible by using X-ray 
micro-tomography to visualize multiphase fluid displacement 
events at the pore-scale within porous media [41,50-54]. The 
images and measurements obtained using this technique have 
contributed towards developing more efficient strategies for 
hydrocarbon recoveries and underground CO, storage [36,40]. 
Recent studies have started to use this technology to investigate 
the pore-scale behavior [33,55] (wettability and trapping) of H3; 
however, these studies have been conducted at ambient conditions 
but not yet at representative subsurface storage conditions. The 
impact of changes in Hz properties between ambient and reservoir 
conditions, such as changes in H2 solubility in water as a function 
of pressure and temperature, requires experimental verification to 
provide an indication of expected fluid behavior during under- 
ground H; storage. 

In this study, we have hence performed a cyclic Hz core-flood 
experiment in a sandstone sample at temperature and pressure 
conditions of 50 °C and 100 bar respectively using a custom- 
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designed flow apparatus, and imaged the core in situ in an X-ray 
micro-CT scanner. The sample was imaged before and after drai- 
nage and imbibition, followed by segmentation and image analysis 
to quantify the phase volumes, H2 saturation, trapping, and contact 
angles. We provide a quantitative analysis of the amount of H2 
trapped in the subsurface during underground Hp storage and 
determine contact angles in situ. Furthermore, we examine the 
potential effect of dissolution of Hz by providing a comparison of 
the residual Hz saturations when using H2-equilibrated brine and 
non-H>-equilibrated brine as the imbibing fluid. 


2. Materials and methods 
2.1. Equipment and materials 


A water-wet Bentheimer sandstone (typical composition of 
Quartz: 91.7%, Feldspar: 4.86%, Clay: 2.68%, Pyrite and Iron 
Hydroxides: 0.17%) [56] sample with a diameter of 6 mm and a 
length of 27 mm was used as the porous medium. Bentheimer 
sandstone is a well-characterized texturally homogeneous rock, 


6.22 mm 
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known to be water-wet with high porosity (20 to 26%) 
[56-58] and permeability (052 to 3.54 Darcy or 


5.13 x 107” to3.49 x 10°"? m?) [56], making it a suitable candi- 
date to initiate the first investigations for pore-scale H2-brine- 
rock interactions at reservoir conditions that can be visualized 
using a micro-CT scanner. The sample’s permeability to water 
was measured separately and found to be 2.9 Darcy 
(2.86 x 10°” m2). 

Before starting the flow experiment, the sample was immersed 
in methanol under a fume-hood for 12 h, and then dried in a vac- 
uum oven under vacuum and at 100 °C for 18 h. A solution of 4 wt% 
potassium iodide (KI) salt (Sigma-Aldrich, UK) in deionized water 
(Thermo Scientific™ Barnstead™ LabTower™ TII Water Purification 
System) was used as the aqueous phase, which provided an effec- 
tive X-ray contrast between brine (wetting phase) and H2 (non- 
wetting phase). High purity (>99.99%) Hz gas (BOC) was used as 
the gas (non-wetting) phase. 

The experiment was performed using the flow apparatus shown 
in Figure S1 (Supplementary Material). It consists of a custom- 
designed high-pressure and high-temperature PEEK core holder 


Fig. 1. 2D and 3D visualization of fluid distributions after drainage with Hz and imbibition with H,-equilibrated brine. (A) 2D vertical cross-section of the filtered image after 
drainage, (B) 2D vertical cross-section of the filtered image after imbibition, (C) 3D visualization of H, after drainage and (D) 3D visualization of Hz after imbibition. The 2D 
images show that after imbibition, the brine ‘B’ replaces H2 ‘H’ in the smaller pores (see for example the pore space inside the yellow marking represented by ‘X’), while H2 
remains in the larger pores. ‘R’ represents the rock grains. The 3D images show that the connected Hz phase after drainage is fragmented into several H2 ganglia after 
imbibition. Red, yellow, and blue colours respectively indicate H3 ganglia occupying pores of sizes that are more than 100 times, between 10 and 100 times, and up to 10 
times larger than the average pore size. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) 


318 


Z. Jangda, H. Menke, A. Busch et al. 


(rs systems), four high-pressure pumps (ISCO, three model 500D, 
and one model 100DX), and a reactor (Parr Instruments Company). 
The sample was wrapped with Teflon tape and aluminum foil and 
placed inside a Viton sleeve connected to the flow lines. The core 
holder was wrapped with a heating tape and securely fitted on 
top of the rotation stage of the X-ray micro-CT scanner (EasyTom 
150, RX Solutions). 

A thermocouple was inserted through one of the ports from the 
base of the core holder and attached outside the Viton sleeve with 
aluminum tape to measure the temperature outside the rock sam- 
ple. To prevent the overheating of the system, another thermocou- 
ple was attached to the heating tape wrapped around the core 
holder. Four high-precision syringe pumps were used to apply con- 
fining pressure, and to inject and receive the fluids. A Hastelloy 
reactor was included in the flow loop to pre-equilibrate the brine 
with Hə for 24 h at the experimental pressure and temperature 
conditions. The reactor was filled with brine up to 95% of its vol- 
ume, with the remaining volume being Hg. The fluids in the reactor 
were then heated to 50 °C and mixed using an entrainment stirrer 
until equilibrated while an ISCO pump provided a constant back 
pressure and a reservoir of H2 gas. The equilibrated H2 and brine 
were then pulled from the reactor’s cell and through the core using 
the receiving pump. This pre-equilibration of fluids was carried out 
to avoid any possible dissolution of H2 in the imbibing brine during 
injection in cycle 1 (discussed in the forthcoming section). 


2.2. Experimental methodology 
The experimental protocol consisted of the following steps: 


i. loaded the sample into the core holder; 

. applied a confining pressure of 20 bar around the sample 
with deionized water; 

iii. flushed the sample with CO, from the base of the sample to 
displace any air inside the sample; 

. Injected 75 pore volumes (PV) (where, 1 PV = 0.19 mL) of 
non-H,-equilibrated brine from the base of the sample at a 
flow rate of 0.5 mL.min“', ensuring 100% brine saturation; 

. switched on the heating to gradually increase the tempera- 
ture of the flow system to 50 °C in one hour; 


(Drainage 
Cimb. eq.brine 
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vi. pressurised the flow system with non-H2-equilibrated brine 
in steps of 10 bar each, while maintaining a 20 bar difference 
between the pore and the confining fluid pressures, to reach 
a pore pressure of 100 bar and a confining pressure of 
120 bar; 

injected 70 PV of H2-equilibrated brine from the reactor to 
completely displace non-H2-equilibrated brine from the 
sample; 

started cycle 1 and injected 15 PV of H2 (non-wetting phase) 
from the reactor to commence the first drainage cycle at a 
flow rate of 0.05 mL.min“! ensuring a capillary number 


(Ca) of 4.2 x 10°°. Here, Ca = oe where ų is the H2 viscosity 
(9.64 x10 Pa- s) [59], V is the velocity of the injected 


fluid (Hz) and y is the interfacial tension (0.0683 Nm~') 
[60] between Hz and water at the experimental conditions. 
In this step, Hz was injected from the top of the sample; 

ix. injected 5 PV of H2-equilibrated brine (wetting phase) from 
the reactor to displace H; at the same flow rate ensuring a Ca 
of 2.3 x 10°°. In this step, brine was injected from the base 
of the sample. The sample was scanned after leaving it over- 
night and it showed complete removal of Hz from the pore 
space, which could be due to diffusion of H2 through the 
Viton sleeve into the confining water along minute and 
microscopically undetectable channels in the wrapping of 
the core sample; 

. started cycle 2 by repeating step viii, i.e., injected 15 PV of Hz 

from the top of the sample to commence the second drai- 

nage cycle with the same Ca of 4.2 x 10°° as the first drai- 
nage cycle. 

injected 5 PV of non-H>-equilibrated brine (wetting phase) 

from the base of the sample using a pump with the same 


Ca of 2.3 x 10°° as the first imbibition cycle. 


vii. 


Vili. 


xi. 


The central vertical section of the sample, which had a length of 
6.7 mm, was X-ray scanned after step (ii) to image the dry sample, 
and immediately after drainage or imbibition in steps (viii) to (xi) 
to visualize the fluids after each displacement process. All scans 
were acquired using an EasyTom 150 X-ray micro-CT machine. 
The scans were acquired with a voxel size of 5 um and with 
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Fig. 2. H) pore occupancy and saturation profiles. (A) Radius of the pores occupied by H, after drainage with Hz and imbibition with H2-equilibrated brine and (B) H3 
saturation profile along the height of the sample after drainage with H, and imbibition with H2-equilibrated brine. 
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1792 projections. The tube voltage was 90 kV and the X-ray tube 
power was 10 W. The duration of each scan was 2 h. 


2.3. Image processing 


The raw data from each scan was reconstructed using the Easy- 
Tom Xact software. The resulting images were then processed and 
visualized using Avizo 2021.2 (ThermoFisher Scientific) as shown 
in Figure S2 (Supplementary Material). A cylindrical mask was 
applied on the reconstructed data to remove the regions outside 
the core sample including the Viton sleeve, aluminum foil and 
Teflon tape wrapped around the sample. A sub-volume of 1250 
slices corresponding to a sample length of 6.22 mm was selected 
for image processing, on which a non-local means filter [61] was 
applied to remove the noise. All wet images (containing Hz and 
brine) were registered to the dry reference image and then resam- 
pled onto the same voxel grid as used in the dry reference image. 
Image segmentation was then performed using a watershed algo- 
rithm [62] based on the grayscale intensity of each voxel. The 
dry sample was segmented into pores and grains which served 


6.22 mm 
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as the mask for the wet segmented images. The porosity value of 
22.8% obtained from the segmented volume is in good agreement 
with reported Bentheimer porosity values between 20 and 26% 
[57,58]. The image segmentation workflow for wet scans (with flu- 
ids) is detailed in Figure S3 (Supplementary Material), and the 
grayscale histograms for the raw and filtered images are shown 
in Figures S4 to S8 (Supplementary Material). 


2.4, Saturation, pore occupancy and contact angle 


Quantitative analysis was then performed on the segmented 
images to calculate Hz saturations, pore occupancy, and contact 
angles. H2 phase volume after each displacement process was used 
to calculate the initial and residual saturations. To obtain the Hz 
pore occupancy after each displacement process, first the pores 
from the dry scan were extracted using a network extraction algo- 
rithm [63]. The algorithm employs the modified maximal ball algo- 
rithm that extracts simplified networks of pores and throats [64], 
by finding constrictions and wider regions in the pore space and 
hence identifies where pores in the image are separated by throat 


Fig. 3. 2D and 3D visualization of fluid distributions after drainage with Hz and imbibition with non-H2-equilibrated brine. (A) 2D cross-section of the filtered image after 
drainage, (B) 2D cross-section of the filtered image after imbibition. (C) 3D visualization of Hz after drainage and (D) 3D visualization of Hz after imbibition. The 2D images 
show that after imbibition, the brine ‘B’ replaces H2 ‘H’ even in some of the pores (see for example the pore space inside the yellow circles represented by ‘Y’) which remained 
saturated with H3 after imbibition with H2-equilibrated brine (compared with Fig. 1B). ‘R’ represents the rock grains. The 3D images show that the connected H3 phase after 
drainage is fragmented into several H, ganglia after imbibition. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of 


this article.) 
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surfaces [65]. The extracted network statistics are provided in 
Table S1 (Supplementary Material). The pores (extracted as 
spheres from the network) were then overlayed on each of the seg- 
mented H2 image, to determine the radii and the number of pores 
filled with H3 after each displacement process. 

To identify the wettability of the system after each imbibition 
step, we measured the H2/brine/rock contact angles on each voxel 
at the three-phase contact line on a sub-volume of 2.5 mm? using 
an automated algorithm [66] with a low standard deviation (© of 8 
and 12 for the two experimental datasets) [66]. This method has 
been used in multiple studies [67-71] to automatically measure 
contact angle values for X-ray tomography images with reasonable 
accuracy. However, using this in situ algorithm at the current 
image resolution, there is a possibility of some error related to 
voxel resolution at the three-phase contact line [33]. 


3. Results and discussion 
3.1. Drainage and imbibition cycle 1 (H2-equilibrated brine) 


The imaging results of Hz drainage and imbibition with pre- 
equilibrated brine are shown in Fig. 1. The filtered 2D vertical 
cross-sections in Fig. 1A and B show the displacement of H2 with 
brine, predominantly in the smaller pores after imbibition. 
Fig. 1C shows a well-connected percolating H, cluster (red) with 
a few isolated Hz ganglia shown in blue and yellow. During imbi- 
bition, Hz fragments into smaller clusters (Fig. 1D). The colours 
in Fig. 1C and D represent the different sizes of the H2 ganglia clas- 


sified based on the average pore size (S) of 1.12 x 10° um?, where 


e Red - H; ganglia occupying pore sizes 100 times larger than S 
e Yellow - H2 ganglia occupying pore sizes 10 to 100 times larger 


than S 
e Blue - H; ganglia occupying pore sizes smaller than 10 times S 


The H saturation calculated after drainage was 36.3%, which 
reduced to 24.8% after imbibition with H2-equilibrated brine, giv- 
ing a recovery of 31.6%. 
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The H2 pore occupancy is shown in Fig. 2A. The shift of H2 pore 
occupancy towards larger pores indicates that the larger pores are 
preferentially occupied by trapped H; after imbibition. The average 
pore radius of the pores occupied by H2 increases from 36 um (o of 
11 um) after drainage to 41 um (o of 12 um) after imbibition. The 
Hə saturation profile along the height of the sample is plotted in 
Fig. 2B. The saturation varies between 24% and 48% during drai- 
nage and 15% and 35% during imbibition, showing some hetero- 
geneity in the sample. Furthermore, the change in saturation 
between drainage and imbibition remains constant over the length 
of the visualized subsection, indicating the absence of capillary end 
effects. The absence of capillary end effects was expected as only a 
6.2 mm long central section of the 27 mm core sample is analysed. 


3.2. Drainage and imbibition cycle 2 (non-H>-equilibrated Brine) 


The pore occupancy for cycle 2 after drainage with H2 and imbi- 
bition with non-H2-equilibrated brine is shown in Fig. 3. The fil- 
tered 2D vertical cross-sections in Fig. 3A and B show that brine 
displaced H; in a larger number of pores compared to cycle 1. Like 
cycle 1, H2 mostly exists as a single cluster after drainage (Fig. 3C) 
which is then fragmented into smaller clusters during imbibition 
(Fig. 3D). The colours in Fig. 3 represent Hz ganglia occupying pores 


of different sizes classified based on an S value of 1.12 x 10° um? 
using the same colour classification scheme as for cycle 1 discussed 
above. 

The H; saturation after drainage was 35.8%. The residual H, sat- 
uration after imbibition with non-H2-equilibrated brine was 20.4%, 
indicating a recovery of 43.1%. However, we posit that it is likely 
that some of the Hə dissolved in the non-H>2-equilibrated brine 
resulting in lower residual saturation in this cycle. Moreover, iso- 
lated clusters of Hz in Fig. 3D appear visually different than those 
in Fig. 1D, which will be compared quantitatively in the forthcom- 
ing sections. 

Fig. 4A shows the pore occupancy of H2 calculated using image 
analysis for both drainage and imbibition in cycle 2. We observe 
the same shift as in cycle 1, that is larger pores remained occupied 
by H2 after imbibition. The average radii of the pores occupied by 
Hə are 36 um (o of 12 um) after drainage, and 43 um (© of 
12 um) after imbibition. These values are similar to the values 
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Fig. 4. Hz pore occupancy and saturation profiles. (A) Radius of the pores occupied by H} after drainage with Hz and imbibition with non-H,-equilibrated brine and (B) H3 
saturation profile along the height of the sample after drainage with H, and imbibition with non-H2-equilibrated brine. 
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obtained after cycle 1; however, the reduction in the number of 
pores occupied by H2 after imbibition is much higher. This differ- 
ence shows that although H; still preferentially occupies the larger 
pores after imbibition, the number of pores occupied by Hz is sig- 
nificantly lower, again indicating possible dissolution of H, in the 
imbibing brine. The H, saturation profile along the height of the 
sample is shown in Fig. 4B. The saturation varies between 24% 
and 47% during drainage and 10% and 28% during imbibition. The 
change in saturation between drainage and imbibition is larger 
than in cycle 1, indicating a greater loss of Hz. However, this 
change is still constant over the length of the sample. 


3.3. Comparison between cycles 


Fig. 5A shows a comparison of the H, saturation along the 
height of the sample for drainage cycles 1 and 2, while Fig. 5B 
shows the Hz pore occupancy. We observe very little difference 
between the two drainage cycles, indicating a repeatability of the 
experimental protocol for drainage. Both, the qualitative and quan- 
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titative analysis show that the two drainage steps are similar com- 
pared to the two imbibition steps. 

The average pore radii and the number of pores occupied by Hz 
after both drainage cycles (6889 after drainage 1 and 6756 after 
drainage 2) are similar. The saturation profiles are also similar 
for both cycles (Fig. 5A). However, we observe significant differ- 
ences between cycle 1 and 2 for imbibition (Fig. 5C). The residual 
Hə saturation in cycle 2 is lower compared to cycle 1, especially 
in the bottom part of the sample. This difference may indicate that 
during imbibition in cycle 2, some H; dissolves into the brine near 
the inlet, lowering the residual H3 saturation in that part of the 
sample. As the non-H2-equilibrated brine travels further into the 
sample, the brine starts to get equilibrated with H2 and dissolves 
less H2 in the upper part of the sample. The saturation profile then 
starts to align with the saturation profile during the imbibition 
step of cycle 1. The pore occupancy profiles for H2 (Fig. 5D) for both 
imbibition steps are similar in terms of the preferred size of pores 
occupied by H3; however, the number of occupied pores (3331) is 
significantly lower in the case of imbibition with non-H>- 
equilibrated brine, compared to 4333 for imbibition with H3- 
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Fig. 5. H2 saturation profiles and pore occupancy distributions. (A) H2 saturation profile along the height of the sample after the two drainage steps, (B) Radius of the pores 
occupied by H; after the two drainage steps, (C) H> saturation profile along the height of the sample after imbibition with 5 PV of H2-equilibrated brine and 5 PV of non-H;- 
equilibrated brine, and (D) Radius of the pores occupied by H; after imbibition with 5 PV of H2-equilibrated brine and 5 PV of non-H2-equilibrated brine. 
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equilibrated brine, which is expected considering the lower resid- 
ual saturation that was observed. 

Image analysis also shows that imbibition with non-H>- 
equilibrated brine (in cycle 2) resulted in brine occupying pores 
which remained filled with Hz in the previous imbibition (in cycle 
1) as shown in Fig. 6. Although the solubility of Hz in water is low 
(~0.001 g.kg~! at the experimental conditions [72]), due to the 
small size of Hz molecules, the mass of Hz occupying the pore 
space is small. Hence, up to 7% of the H2 present in the pores during 
the injection of 5 PV of non-H2-equilibrated brine could potentially 
dissolve. This dissolution could explain the lower residual satura- 
tion after imbibition with non-H2-equilibrated brine. 


3.4. Comparison of fluid saturations 


Quantification of fluid saturations using in situ X-ray flow visu- 
alization experiments has been frequently reported in the litera- 
ture for CO, Nz, and oil-brine systems. Some of the earlier 
reported studies for CO 2-brine systems on sandstone samples 
{37,40,73-75] are plotted in Fig. 7 and compared with our results. 


Fig. 6. 2D X-ray CT image after imbibition with non-H2-equilibrated brine. The 
bright spots (B*) highlight pore spaces that were filled with brine after imbibition 
with non-H,-equilibrated brine in cycle 2, which were occupied by H, after 
imbibition with H2-equilibrated brine in cycle 1. 
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Fig. 7. Initial and residual H2 saturations for the two imbibition steps in this study 
compared with literature data for Hə and CO, obtained from X-ray flow visualiza- 
tion experiments. 
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Four of these earlier studies [37,73-75] are at experimental condi- 
tions similar to those used in this work and the recovery trends are 
qualitatively similar. Interestingly, the initial and residual satura- 
tion values of CO; closest to our results are for the reported values 
at ambient conditions (with CO, as the gaseous phase) [40]. This 
could be due to the relatively smaller difference in the density of 
Hə at the temperature and pressure used in our experiment and 
CO, at ambient conditions, compared to CO, at high temperatures 
and pressures. 

Recent experimental measurements for Hz saturation in a H2- 
brine system [55] at ambient conditions are also plotted in Fig. 7. 
In this experiment, an initial Hz saturation of 64% and a residual 
Hə saturation of 41% after 5 PV of brine injection was reported 
for a Gosford sandstone [55]. Although the recovery trend of our 
results broadly agrees with this work, the reported initial and 
residual Hz saturation values are significantly different from our 
results. This difference could be due to the higher temperature 
and pressure conditions and a larger sub-volume analysed in our 
work, which provides a more realistic scenario for subsurface 
storage. 

Additionally, the decrease in H2 residual saturation due to dis- 
solution of Hz in the non-H2-equilibrated brine is likely to be an 
important parameter for reservoir-scale modelling, especially dur- 
ing cyclic injection of H2. The initial recovery from the first injec- 
tion into the reservoir may be low with subsequent cycles 
producing more H as the reservoir brine is slowly equilibrated 
and loss of Hz gas due to dissolution and trapping decreases. The 
impact of dissolution and other important factors like microbial 
activity and diffusion on the fluid saturations during underground 
Hə storage warrants further investigations. 


3.5. Contact angle 


Contact angles were determined for a 2.5 mm? sub-volume for 
each imbibition step. The results (Fig. 8) show water-wet condi- 
tions with an average contact angle (measured through brine) of 
54° after imbibition with H2-equilibrated brine (cycle 1) and 53° 
after imbibition with non-H2-equilibrated brine (cycle 2). These 
contact angle values agree with the measurements performed 
using multiple methods on a H2-brine system in a Bentheimer 
sandstone [33], where average contact angles were reported to 
be in the range 39.8° to 59.8°. The mean contact angle (59.8°), 
which was obtained via the same 3D local method [66] used in 
our study, agrees well with our results. Other studies [24,48,49] 
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Fig. 8. Contact angle distributions after imbibition with H -equilibrated brine 
(cycle 1) and non-H>-equilibrated brine (cycle 2). 
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have found sandstone samples to be strongly water-wet for H2- 
brine systems. They have reported contact angles between 6° and 
43° regardless of the tested parameters for temperatures between 
20 °C and 50 °C and pressures between 20 bar and 100 bar. 


4. Conclusion 


In this study we have imaged the spatial distribution and trap- 
ping of H, after the injection of both H2-equilibrated and non-H>- 
equilibrated brine in a Bentheimer sandstone at subsurface pres- 
sure and temperature conditions. The analysis shows that the rock 
sample is water-wet at subsurface conditions with an average con- 
tact angle of 54 and 53 obtained after H2-equilibrated and non-H>- 
equilibrated brine injection, respectively. These observations are in 
agreement with some of the recent experimental studies [33,35]. 
Other studies [24,48,49] have also reported strongly water-wet 
conditions for H2-brine systems. These water-wet conditions can 
enhance snap-off processes resulting in more trapping of H2 [76]. 
The robust determination of wettability is hence key to under- 
standing and quantifying trapping and recovery. More in situ 
experimental datasets are needed for different pressure and tem- 
perature conditions as well as rock types and flow rates to under- 
stand how the pore-scale displacement processes impact large- 
scale H, storage. 

Our analysis presents new insights into the possible dissolution 
of H3 into brine at subsurface conditions. The results show a 
decrease in residual H2 saturation and a decrease in H2 pore occu- 
pancy when the sample was imbibed with non-H,-equilibrated 
brine, indicating that dissolution of Hz may be an important con- 
sideration for the accurate prediction of Hz recovery. Since in situ 
experimental work on H; at high temperature and pressure condi- 
tions is still scarce, this effect was not reported in any previous 
work but can be a critical aspect to consider for future research. 

We have found that the trend in recovery factor for our exper- 
iment is consistent with previously reported observations [37,73- 
75] for supercritical CO2-brine experiments at similar pressure and 
temperature conditions. However, the residual saturation values 
(corresponding to the initial saturation values) closest to our 
results are for CO2-brine experiment at ambient conditions [40]. 

In summary, even though porous geological formations con- 
taining depleted gas reservoirs or saline aquifers can provide large 
storage volumes for Hz, pore-scale interactions at reservoir condi- 
tions need to be further investigated to provide the key physical 
insights if and why these porous reservoirs are suitable for long- 
term and cyclic underground storage. A key technique is this 
regard is utilizing fast synchrotron X-ray imaging [52,53,67,77]| 
to obtain time-resolved 3D images for visualization and quantifica- 
tion of the dynamics of fluid displacement during H3 injection and 
production. 
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